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SUMMARY 

This work has Investigated the range of accuracies achieved In the gear 
tooth temperature using an automatic finite element mesh generator. Gear web 
contribution to the gear cooling process Is studied by Introducing a varying 
size hole at the center of the gear because of the versatility of program TARG 
In allowing different heat transfer coefficients In different areas of the 
gear tooth. A study Is carried out to evaluate the contribution of the loaded 
and unloaded faces as well as the top and bottom lands. 

A general purpose two-dimensional finite element preprocessor ATOGEN has 
been developed for automatic generation of a finite element mesh over a pie- 
shaped sector of a gear. The program Is used for facilitating the Input to an 
upgraded version of a previously developed program for the thermal analysis of 
running gears (TARG). The latter program determined the steady state tempera- 
ture distribution throughout the specified gear. The automatic mesh generator 
program Includes a band width minimization routine for reducing computer cost. 


INTRODUCTION 

The problem of predicting gear tooth temperature by means of a finite 
element technique has been recently treated by several authors (refs. 1 to 3). 
The finite element mesh used In these references, though applicable for any 
gear tooth proportions. Including nonstandard teeth, was limited to a fixed 
number of elements along the tooth profile. Eight (8) divisions along the 
tooth appeared to be adequate at the time for predicting gear tooth temperature 
with reasonable accuracy. Because of the sharp temperature gradients observed 
In the results. It became necessary to develop a finite element preprocessor 
capable of descretlzlng the tooth Into any required mesh. In addition to the 



customary element connections and nodal coordinate lists that an all purpose 
preprocessor would yield, boundary nodes on loaded and unloaded sides as well 
as top and bottom lands would be accordingly Identified. 

Parameters which are pertinent to the gear geometry such as all nodes on 
the Involute profiles, radii of curvatures at the same, and adjacent nodes to 
the pie-shaped boundary constitute necessary output for executing the thermal 
analysis code. The Program ATOGEN has been developed for this purpose. The 
program optionally punches an output file for subsequent Input to the Thermal 
Analysis of Running Gears program "TARG" which ultimately determines the steady 
state temperature distribution throughout the specified gear. Parametric 
studies are carried out to evaluate the effects on tooth temperatures of the 
mesh size (number of segments along tooth profile), the significance of gear 
blank size, and the Individual contributions of the tooth boundaries to the 
convection cooling process. 


DESCRIPTION OF PROGRAM ATOGEN 

The program Is partly based on methods developed In references 4 and 5. 
The pie-segment of a given gear Is Initially divided Into a crude mesh made up 
of quadrilateral (four-sided) regions as shown In figure 1. The sides of each 
Individual region are generally curvilinear 1 and a triangular region can be 
viewed as a quadrilateral by dividing one of the sides Into two segments. The 
corners of the four sided regions In the tooth portion of the model Identify 
significant points on the tooth profile. Starting from the outermost edge of 
the tooth these are the points on the top land, the pitch points, the lowest 
points of tooth contact and the points on the bottom land. Intermediate, or 
mid-side nodes are added on each side of the quadrilateral for achieving the 
closest approximation of the gear Involute profile. 

The coordinates of the grid points shown In figure 1 are computed in the 
program for any gear geometry having a modified Involute tooth profile with a 
given number of teeth, diametral pitch, pressure angle, profile shift coeffi- 
cient, whole depth, and topping. The x and y coordinates of grid points 
22 through 28 are given by 

x. = r. sin e. 

1 1 1 1 = 22, 23 28 (1) 

y 1 = r i cos e 1 

where 

r^ = radius at grid point 1 
= angle from the y-axls to the radius at 1 
Using Involutometry, the angle can be written as 


^or simplicity, In the present version of the program, the sides are 
assumed as quadratic curves. 
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( 2 ) 


= Inv (D) + ]/4Pd + Cp tanD/Pd 2 - Inv 1 O^/r ^) 2 “ 1 

where 

$ = pressure angle 
P(j = diametral pitch 
Cp = pressure shift coefficient 
rb = base radius 

A conformal mapping of the general quadrilateral In the x - y plane 
Into a square In the 5 - n plane Is accomplished using transformation 


N 


= £ x^ , y = 53 ^ 


(3) 


1=1 


1=1 


where. 


N r = number of nodes around the quadrilateral 

For an eight-node quadrilateral, the shape functions are given as 
(see ref. 6 ) 


N n =}(!- 5)0 - n)(5 ♦ n ♦ 1 ) 


N 2 = 2 (1 “ 5)(1 ' n) 


N 3 = 4 0 + OO - n)(5 - n - 1) 
N 4 = \ (1 + 5)0 - n 2 ) 


N 5 = j 0 + 5)0 + M(5 + n -1 ) 


N 6 = \ O - 5 2 )0 + n) 


“ 4 O - 5)0 + n)(n -5-1) 


N 8 = 0 - 5)0 - n ) 


\ 


J 


(4) 


The pair (x, y) represents the coordinates of node 1 on the quadrilat- 
eral. The transformation Is shown schematically In figure 2. 
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Once the transformation Is accomplished, the square region In the i - n 
space Is divided equally Into a number of rows (parallel to the ij-axls and 
a number of columns. To avoid duplication of nodes, an equlvalenclng procedure 
Is undertaken by assigning a connectivity matrix JT(I, J) to each region I 
with sides J as shown In figure 3. Nodes on adjoining sides are checked for 
duplicity via JT(I, J) and would acquire their original numbers If duplication 
Is encountered. The node numbering follows a sequence that starts at 
5 = -l , n = +1 and proceeds from left to right (ref. 4). 

A boundary Identification matrix JB(I, J) Is assigned for each side J 
of a region I to associate the element boundaries with the original segments 
of the overall boundary. The purpose of this Identification matrix Is to dis- 
tinguish areas along the boundary according to the convection heat transfer 
coefficient to be assumed along the exposed portion of the model boundary. 
Specifically two different heat transfer coefficients are assumed along the 
loaded side of the tooth one above the pitch points and one below It, one co- 
efficient for the top land, and another one for the unloaded side of the tooth. 
These heat transfer coefficients are usually derived experimentally or from 
empirical formulas, (see ref. 1). Figure 4 shows the boundary Identification 
numbers (chosen arbitrarily) for region sides that lie on the boundary. Sides 
which are not part of a boundary assume a zero JB value. 

Each region Is then divided Into triangular elements by joining the 
shortest diagonal In the elementary quadrilateral. In addition to providing 
the capability of changing element size within a given region by assigning the 
number of rows and columns. It Is possible to further refine the mesh by 
changing the relative locations of regional nodes on the line of symmetry. 

This Is done because solution accuracy of stretched triangular elements Is 
usually poor and one should strive to approach an equilateral triangle. Once 
the entire region has been descretlzed and node numbers, Individually assigned, 
a new node renumbering scheme Is used to minimize the conductivity matrix 
bandwidth. 

The routine scans for all nodes In the proximity of a given radius, pro- 
ceeding from left to right. Using this approach It Is found that a substantial 
reduction of the bandwidth Is achieved leading to significant savings In com- 
puter time. 


INFLUENCE OF MESH SIZE ON ACCURACY OF DETERMINED TOOTH TEMPERATURES 

The test case evaluated here Is the same one used and experimentally ver- 
ified In reference 1. It consists of two Identical eight pitch standard spur 
gears with 28 teeth each and 20° pressure angle. The face width of each gear 
Is 0.25 In., and the root radius Is 1.575 In. 

Several meshes were generated using program ATOGEN for subsequent process- 
ing through the thermal analysis program TARG. Computer generated plots of 
the meshes In ascending order of number of segments (or levels) across the 
tooth are shown In figures 5 and 6 . The coarsest mesh Is similar to the model 
used In references 2 and 3 and consists of eight segments across the tooth 
profile. The finest mesh shown (enlarged for the tooth portion of the model 
for clarity) has 48 segments. The temperature contours (Isotherms) are shown 
In figure 6 for two extreme meshes of 8 and 48 segments. Although little 
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difference Is apparent In the contour shapes of the two meshes, the average 
temperature could be off by as much as 12 percent when a coarse mesh Is used. 
The trend Is Illustrated In a plot, (fig. 7), of temperature versus number 
segments. As expected, there Is a number of segments beyond which no Improve- 
ment In accuracy Is discernible. Three temperature norms are chosen for this 
Illustration: namely the maximum, average midplane, and average face tempera- 
tures. From the graph It Is concluded that 24 segments will yield an accuracy 
better than 1.0 percent and with 18 segments the error In maximum temperature 
should not exceed 2.5 percent. 


GEAR BLANK CONTRIBUTION TO THE COOLING PROCESS 

Because of the option available In program ATOGEN of modelling a gear 
with an annulus, gear models are generated with various size holes to evaluate 
gear blank contribution to the cooling process. In addition to the baseline 
case with no holes, four other configurations are examined. The finite element 
models for the latter are shown In figure 8. The number of tooth segments 
used In all models Is 18. A plot of the resulting temperatures versus hole 
size Is shown In figure 9. The graph shows that the gear blank has very little 
Influence on tooth temperatures and the heat transfer process Is limited to 
the teeth regions. The steady state heat transfer problem In gearing can 
therefore be analyzed by considering only the tooth portion of the gear with 
little or no loss In accuracy. 


BOUNDARY COOLING CONTRIBUTIONS 

Included In the printout of a program TARG run Is a heat balance check 
which assigns heat loss contributions of the different portions of the model. 
The convective heat transfer coefficients used In the present case are the 
same as those noted In reference 1 with nomenclature Illustrated In figure 
10. The percentage heat loss attributed to the lubricated flank of the tooth 
Is 24.3 percent, that to the unlubricated flank Is 43.1 percent, that to the 
gear transverse sides Is 26.0 percent, and 6.6 percent to the top land of the 
tooth. Unfortunately these allocations are heavily dependent on the numerical 
value assumed for the heat transfer coefficients and would require experimental 
verification. 


SUMMARY OF RESULTS 

A brief description Is given of the finite element preprocessor program 
ATOGEN developed for the automatic generation of finite element meshes of 
gears. The program Is designed to produce Input data necessary for running 
the steady state heat transfer analysis program TARG and Includes In Its code 
gears with modified Involute profiles. In addition to the punch files, op- 
tional plots of the generated mesh can be produced by ATOGEN. 

The advantages gained from having a finite element preprocessor are quite 
numerous. Some of these are self evident such as the time savings over manual 
preparation of element, nodal, and boundary data. Computer cost can be sub- 
stantially reduced by judiciously assigning a finer mesh to areas of higher 
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temperature gradient and using a coarse mesh elsewhere. Cost Is also minimized 
through the use of a matrix bandwidth minimi zer. These advantages were uti- 
lized In the present Investigation with the following results obtained: 

1. It was found that substantial errors in the predicted tooth tempera- 
tures could result from using a coarse mesh, such as the one used In reference 
1. For the example shown the error In maximum temperature approached 12 per- 
cent. This error could be held to less than about 1 percent by using a mesh 
with 24 or more segments across the tooth. 

2. The amount of heat dissipation through the gear blank was fond to be 
negligible for the cases studied. The error Incurred In the temperatures by 
considering only the tooth portion of the gear Is less than 3 percent. 
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Figure 2. - Transformation coordinates. 
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Figure 3. - Regional connectivity matrix layout. 










Figure 5. - Computer generated mesh. Enlarged plot 
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